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Abstract 

The plasmon excitations in proposed single- and double-component helical liquid (HL) mod- 
els are investigated within the random-phase approximation, by calculating the density-density, 
spin-density and spin-spin waves. The effect due to broken time-reversal symmetry on intraband- 
plasmon dispersion relation in the single-component HL system is analyzed and compared to 
those of well-known cases, such as conventional quasi-one-dimensional electron gases and armchair 
graphene nanoribbons. The equivalence between the density-density wave in the single-component 
HL to the coupled spin-density and density-density waves in the double-component HL is shown 
here and explained, in addition to the difference between intraband and interband-plasmon excita- 
tions in these two systems. Since the two-component HL can physically be thought of as a Kramers 
pair in two-dimensional topological insulators, our proposed single-component HL model with bro- 
ken time- reversal symmetry, which is an artificial construct, can be viewed as an "effective" model 
in this sense and its prediction may be verified in realistic systems in future experiments. 
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I. INTRODUCTION 



Topological insulators (TIs) are found to be a new class of materials which possess in- 
sulating (large bandgap) states in the bulk but conducting edge states on their surfaces. 
Only states along edges of TIs have non-fiattened dispersion (extended states) within a 
bulk bandgap, thereby allowing for charge and spin currents under zero bias [1]. The bulk 
bandgap of a TI is usually large enough to be comparable to room temperature, which 
makes TIs thermally stable and a good candidate for high-power electrical and optical ap- 
plications. One example of such TIs is the so-called three-dimensional (3D) TI, e.g., Bi2Te3 
crystals. These 3DTIs are 3D band insulators possessing two-dimensional (2D) conducting 
surface states. Additionally, there exist 2DTIs, in which the conducting states are localized 
close to all edges of a slab in real space and display a well-defined Dirac cone in momentum 
space. In this paper, we confine our attention to another predicted type of TI which has 
already been demonstrated in an inverted HgTe/CdTe (also called type-Ill) quantum well 
when the well thickness exceeds 6.3 nm. The type-Ill quantum well is formed by sandwich- 
ing a thin narrow-bandgap HgTe layer between two thick wide-bandgap CdTe layers in the 
growth {z) direction. Bulk HgTe/CdTe is a semimetal with zero bandgap. The quantum- 
size effect from a quantum well introduces a finite but very small bandgap to a HgTe/CdTe 
layer. As expected, the induced bandgap in a HgTe/CdTe layer decreases with the layer 
thickness. For an inverted HgTe/CdTe quantum well, the usual lower p-type Fg-band moves 
above the s-type Fg-band around the center of the Brillouin zone to create an anti-crossing 
bulk bandgap as well as a negative effective mass for conduction electrons at the same time. 
This HgTe/CdTe-based TI is one type of 2DTI. A semi-infinite quantum well extends in- 
finitely in the x direction but is still confined within a half-plane {y < 0) in the y direction. 
Consequently, there always exists one conducting channel (helical edge state) for spin-up 
or spin-down electrons along the edge {y = 0) of a half-xy plane (with a finite thickness in 
the z direction). This means that we may easily assign a spin component to a current if 
we know its flowing direction. The switching from a left to a right circular-polarization of 
incident light is expected to change the flow direction of a photocurrent [2] . This helical state 
spans over a quantum well in the z direction and forms a quasi-one-dimensional electron 
gas (quasi-lDEG). The energy dispersion of the helical state lies within the anti-crossing 
bandgap region of an inverted HgTe/CdTe quantum well. The theoretical description of this 
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type of topological states is given by the model of Bernevig, Hugues and Zhang (BHZ) [3l H] . 
The general classification of TI is present in Refs. [T^ fIU\ . 

The spin dynamics of TIs have received a great deal of attention, including a topolog- 
ical quantum phase transition in a tunable spin-orbit system [S], spin-polarized electrical 
current |6] and a photocurrent induced by circularly-polarized incident light [2j. However, 
unique properties of charge dynamics in the same systems are much less known. The goal of 
this paper is to investigate the charge density/spin dynamics of collective excitations of those 
edge-bound electrons. To study the charge dynamics of undamped collective excitations of 
electrons having a given helicity, we employ a 2DTI model system [HI H] . Those readers in- 
terested in the collective response of 3DTIs may start by looking up Refs. [HI E] • 3DTIs are 
beyond the scope of the present paper. For the 2DTI system of class All (QSH edge), the 
BHZ model predicts the existence of both insulating bulk and conducting (1, 0)-edge states 
in an inverted HgTe/CdTe quantum well, similar (but not identical) to those in a metallic 
armchair graphene nanoribbon (ANR). In a graphene ANR, the plasmon excitations are 
excited by interband transitions only [TT]. Unlike those in a 2DTI system, electronic states 
in a graphene ANR are not localized around the two ribbon edges. Moreover, both heli- 
cal branches contribute to plasmon excitations. This makes the plasmon dispersion in a 
graphene ANR almost identical to that of a conventional IDEG. 

In this paper we investigated the effect of wave function localization in a semi-infinite 
2DTI system on the collective excitation of electrons. Our calculated plasmon dispersion is 
compared with those in a conventional IDEG [12] and in a metallic graphene ANR. Following 
the work by Zhang et al. [TH], we adopt the terminology that an n-component helical liquid 
(HL) contains n-time reversal pairs of fermions. These n-component HL states are localized 
on different edges of a quantum well. For a semi-infinite quantum well, on the other hand, 
we may only need to consider a single mode of the pair to calculate the Coulomb excitation 
of electrons. Since throughout the paper we do restrain ourselves to the semi-infinite case 
we shall refer to such single edge mode as single component HL. Since its counterpart is 
ignored the time reversal symmetry is broken. That is we have a quantum Hall state. Recall 
that the quantized charge Hall (QH) conductivity is attributed to the transport via a single 
chiral edge mode. 

On the other hand when we do include finite spin into consideration we would have 
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quantum spin Hall (QSH) effect which does not require time reversal breaking. A pair of 
states per edge appear. Here we refer to such case as two component HL thus effectively 
considering Z2-trivial TI[3J. The contributions from a single and two component HL to 
the plasmon excitation are explored. Our calculations indicate that ANR and the semi- 
infinite 2DTI system can be considered to be two and single component HL, respectively, in 
calculating the charge dynamics of collective excitations. 

The rest of the paper is organized as follows. In Sec. [11} we calculated the edge-localized 
helical states of electrons and their energy dispersion in a semi-infinite inverted HgTe/CdTe 



quantum well. In Sec. |III[ the dispersion of edge-plasmon excitations of single component 
HL with broken time-reversal symmetry is calculated. That type of density-density response 
is compared with that of a conventional one-dimensional electron gas, a metallic armchair 
graphene nanoribbon and finally with the response of the Kramers pair of two-component 



HL. In Sec. IV we argue that the density-density plasmon excitations of single component 
HL are formally equivalent to the interference pattern between spin-density and density- 
density plasmons in two component HL Finally, the conclusions of the paper are given 
briefly in SeclVl 



II. FORMALISM OF 2DTI IN HgTe/CdTe QUANTUM WELL 

To study the collective electronic excitations, let us start with the BHZ model for the 
electron band structure near the center Tk = (0, 0) of the first Brillouin zone. We assume 
that the quantum well is infinite along the x direction and finite or semi-infinite along the 
y direction. The width of the well is given implicitly within the parameters [A, B and A 
in Eq. ([T])] of the model Hamiltonian. The "ansatz" wave function is taken to be exp(iA;a;) 
in the x direction. Along the y direction, we discretize the spatial position as ?/ = ja > 
with j = 1, 2, . . . , being a positive integer and a being the lattice constant, where A^ is 
the total number of sites assumed for numerical simulations in the y direction, and is taken 
sufficiently large to ensure that two edge topological states do not overlap each other. The 
wave number k in the x direction is given in units of ks = it /a, leading to discrete spatial 
positions along the x axis as x = ia with i = —N/2, . . . , —2, —1, 0, 1, 2, ... , A^/2 being an 
integer. 

According to Ref.[3l the Hamiltonian describing electronic states in a HgTe/CdTe (type- 
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FIG. 1: (Color online) Full band structure e{k) of a 2DTI with A = 25, ^ = B = 1 and = 30. 
The blue curves correspond to the bulk modes while the two green curves represent the topological 
surface states. Each state is doubly degenerate with respect to the two eigenvalues of the Ux matrix. 

Ill) quantum well can be written as a block-tridiagonal matrix 
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e = A sin(A;) Fi + [A - 45 + 2B cos{k)] F5 



t 



Here, the elements of the Clifford algebra are expressed in terms of Fi = 



cr^, F9 



—ay ® ctq and F5 = cr^ g) ctq with (Xj denoting the Pauli matrices. Additionally, A, 5, A in 
Eq. ([1]), which are scaled by fiksVp with vp ~ c/1000 being the Fermi velocity expressed in 
terms of the speed of light c, are the material parameters dependent on the quantum well 
width. The form of the Hamiltonian in Eq. ([T]) implies that the associated wave function 
vanishes at edges. 

The calculated energy dispersion corresponding to Eq. Q is presented in Fig.[T| For 
< A/S < 4, two localized states at the boundaries are obtained, which may be expressed 
in terms of a linear combination of the ansatz wave functions [141 
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where the hehcal states are chosen to be the eigenstates defined by Tiji) = ± |±). The 
analytic solutions of Eq. ^ may be obtained by substituting the above ansatz wave functions 
into Eq. ([!]) and considering the identity [ir5r2, Fi] = 0. This yields 



i^fik) = (c> Jpi(fc)]^- + cl,[p,ik)y) 1+) (3) 
^fik) = (c< Jpi(-fc)]^+i-^- + c<2[p2(-fc)]^+^-0 I-) 

The wave function corresponding to |+) in Eq. ^ is related to the energy dispersion 
Ej^{k) = AsiiD.{k), while that corresponding to |— ) is associated with E^{k) = —Asin{k), 
as shown by two green curves in Fig.[TJ The deviation from the sine function becomes 
significant once E±{k) merges with the bulk modes (shown as blue curves). Parameters 
Pi{k) and P2{k) in Eq. ^ are defined by 



PiAk) = , (4) 

n{k) = A - 2B [2 - cos{k)] . (5) 

For A = B = 1 and A; = 0, the variation of these parameters with respect to A/S is displayed 
in Fig.[2j 

Since our interest is limited to calculating the plasmon dispersion in the long-wavelength 
limit, we would only consider the case with A; ^ 1. This leads to the approximate expression 
E±{k) ^ ±Ak, and pi^2{k) ~ Pi, 2(0) become fc-independent at the same time. For A = 25, 
we find pi^2 = ±Vv^- 

The coefficients c±^a with a = 1, 2 in Eq. (|3| may be determined from the boundary 
conditions as well as the wave function normalization condition. By choosing ^ 1, 
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FIG. 2: (Color online) Based on Eq. ([s]) the real (solid curves) and imaginary (dashed curves) 
parts of pi (green) and p2 (blue) at A; = as a function of A/S for A = B = 1. 

the condition for a finite-valued wave functions requires „ = c^j: ^ = 0. Moreover, the 
vanishing boundary conditions lead to + 2 = ^ind ^ + cf 2 = 0- After normalizing 
the wave functions in Eq. (|3]), we find 1 = c"^ 1 = 2/-\/3. Therefore, the two components 
of the wave function in Eq. ([2]) do not overlap and are localized on opposite boundaries. 
In addition, the eigenvalues of the spin state of ax are related to the sign of the group 
velocity, i.e., ipj corresponds to E+{k) while is related to E-{k). This pair of fermions 
constitutes a 1-component HL, and is connected by the time-reversal symmetry. For a semi- 
infinite quantum well with — > 00 and pi^2 = ±Vv^, on the other hand, we are left with 
a single helical state localized at the j = boundary 



E+{k) ^ Ak . 

If the material parameters |pi,2| > 1 are chosen, one would retain the left mover i/j^ pro- 
portional to |— ) as the proper solution. For a 2-component HL (Kramers pair), we must 
consider two helical states on each edge. Those two states are related in Eq. ^ by the spin 
change |+) ^ |— ) and the time reversal k ^ —k as well. Therefore, the second part of the 



Kramers pair is: 



E^{k) ^ -Ak . 



III. PLASMONS IN HL COMPARED WITH IDEG 



Based on the calculated full band structure from the BHZ model, we will further study 
the electron screening dynamics from the dielectric function of a 2DTI system. The edge 
states in such a system appear as the Kramers pair (Eq. ([6])) with the electron spin attached 
to its momentum. On the level of density-density response it is imposable to break up the 
Kramers pair into its components, thus observing the response of the two component HL. 
However in the next section we shall show that the density-density response of the single 
component HL formally correspond to the interference pattern of the density-density and 
spin-density waves of the two component HL. 

We shall formally start with the density-density response of the single component HL. 
Specifically, we will consider a semi-infinite type-III quantum well, in which only one helical 
state can occur and is localized around the edge {y = 0) of the system. As a result, there exist 
only intraband transition in our system. For the wave function given by Eq. ([6]), the plasmon 
excitation dispersion u}p{q) within the random-phase approximation (RPA) is determined by 
the zero of the following dielectric function 



e{q,u) = l-Viq)U+,+ {q,u + tO+) (8) 

= l--j^^^ + tqViq)5{hw-Aq) , 
IT nio — Aq 

where the noninteracting polarization function at zero temperature is given by 



, , 1 9(k + q) -e(k) 

_ 2 q 

2n hu - Aq + iO+ 

Here, 6{x) is the unit step function and Ep = is assumed for the Fermi energy. The 
perfector of two in Eq. (|9| accounts for the double degeneracy of the bands. Very importantly, 
the time-reversal symmetry is broken in Eq. ^ for the response function, i.e., H+ ^(g, tu + 
iO^) 7^ H^_,_(g, —u + iO^). The Coulomb matrix element introduced in Eq. (j8| is given 
byHa 
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FIG. 3: (Color online) The q dependence (in units of the reciprocal inter-atomic spacing a) 
of the scaled Coulomb matrix element y(g)/(4e^/esa) in Eq. (10) evaluated between its upper 



(3/40) KQ{q) (green) and lower bound —(3/40) ln(g) (black) for various values of the ratio A/B, 
where A = B = 1. The result for the approximate expression in Eq. (11) is given by red dots in 
the figure. 



(10) 



where Kq{x) is the modified Bessel function of the second kind, e^j = eoCb and is the 
dielectric constant of the host material. If it is not explicitly stated, our parameters of 
choice are A = B = 1 and A = 2B in this paper. For this parameter choice, the summation 



in Eq. (10) can be carried out explicitly to give 



V{q) 



2e' 



12cos(gt) dt 



esa J 5VlT¥ [41 - 9cos(2gt)l 



e,a 20 



(11) 



Here, the validity of the last approximation is demonstrated in Fig.|3| and the perfector 
can be calculated, by using the fine-structure constant and eb = 7 for HgTe, as e^ks/^sA. = 
Att X 0.001{hkBVF/A){c/vF) ~ 47r. 



The imaginary part of the dielectric function in Eq. ([8j) yields the boundary hup^h = Aq 
for the particle-hole excitation region. Collapse of the particle-hole region into the single line 
is characteristic to the linear dispersion [11] . Zeros of the real part give rise to the plasmon 




FIG. 4: (Color online) Scaled plasmon excitation energy ujp{q)/ujQ of 2DTI (blue solid curve) 
given by Eq. (|12|) with A = B = I. The lower blue solid curve corresponds to A = 2B, while the 



upper blue solid curve is for A = B oi A = 3B. The result of Eq. (13) is given by the dashed blue 



curve in the long- wavelength limit. The conventional IDEG plasmon dispersion is indicated by a 
red dashed curve. The particle-hole excitation region is represented by a thin black straight line. 

dispersion relation: 



hup{q) = Aq 



(12) 



In the long-wavelength limit with qa ^ 1, we obtain Re[n+_+(g, u + ^0"*")] = q/nhuj + 0{q^) 



as a leading-order result. In conjunction with Eq. (11), the above equation leads to the 



explicit plasmon dispersion relation in the long-wavelength limit: 



Up{q) = -u,q\Ti{q) + 0{q'') , (13) 
nuin 
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Results from our calculations based on Eqs. (11), (12) and (13) are presented in Fig. 4 It is 



found that A = B = 1 and A = 2B gives us the lower boundary of the plasmon excitation 
energy. This implies that the variations of A, B and A within the regime in which the 
topological edge state exists results in an increase in the plasmon energy. 

Most inelastic scattering experiments measure the dynamic structure factor ~ 
Ini[e~^(g, uj)] or the inverse dielectric function. We display Im[e~^(g, u)] in Fig.js] as a 
function of u for chosen values of q. Clearly, the spectrum of Iin[e~^(g, u)] is domi- 
nated by the plasmon resonance. The particle-hole excitation is not pronounced in this 
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FIG. 5: Spectral function Im[e ^(g, uj)] of single component HL in 2DTI as a function of lj/ijJq 
with and chosen q. Here, the curves with different values for q are offset vertically for clarity. 

figure. Around the plasmon resonances, we can employ the plasmon-pole approximation 
Ini[e^^(g, uj)] ~ Pq5{uj — Up), where the plasmon weight is defined by 



TT 

^'^ = [dRee{q, ^^^^ 



7i^hw^qln'^{q)_ ^^^^ 



Viq) 

Although varying the parameters of a 2DTI system from those producing the minimal UJp{q) 
can increase the plasmon energy, it reduces the plasmon weight. 

Now, let us compare our result with several known cases of IDEG. For a conventional 
semiconductor quantum wire with a parabolic energy dispersion -Eo(^) = h'^k'^/2m* for 
conduction electrons, the intraband plasmon dispersion in the long-wavelength limit can be 
written as [T2l [T5| [16] 



Up{q) = uo gV-ln(g) + 0{q') , (16) 



where uid and m* denote the electron linear density and effective mass, respectively. For 
IDEG, the wave vector q is scaled with the characteristic size W of the nanowire, i.e., 
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FIG. 6: (Color online) Plasmon weight (3q for Tl-based (blue solid curve) and conventional (dashed 
red curve) IDEGs, where ujq is chosen as the same for both cases. 

q — )■ qW . The linear scaling of the plasmon frequency with the square-root of electron 
density as well as the high sensitivity to the wire characteristic size ^ 



In(giy) are the 

unique properties of IDEG in conventional semi-conducting nanowires. These properties are 



in sharp contrast with our result in Eq. ( 16 ) for localized quasi-lDEG in a 2DTI system. 



For the conventional IDEG, the plasmon weight is given by 



/3a 



q [-ln(g)] 



3/2 



2e^ 



[KM + 1.972] 



(17) 



which has different power dependence for the term ~ ln(g). Comparison of the plasmon 
weights are shown in Fig.[6] for both the Tl-based and conventional IDEG. From Fig.|6| 
we find that the plasmon weight in 2DTI is an order of magnitude larger than that of the 
conventional IDEG. Consequently, a much more pronounced dynamical structure factor is 
expected for 2DTI [compare Fig. [s] here with Fig. 7 of Ref. [I2]]. We also note that the 
particle-hole excitation spectrum exists in a wide region kpq/m* — h^^Eo^q) < Up^h < 
kpq/Tn* + h~^EQ{q) for a conventional IDEG rather than a narrow line in 2DTI. 

One of the striking features of IDEG is that the RPA becomes exact for small value 
vectors q. This implies that the quasi-particles obtained from the exactly solvable linearized 
Tomonaga-Liittinger model agree with the calculated plasmon dispersion relation in the 
RPA |T7|. In the case of a single-component HL for a 2DTI system, only one electron branch 
in the Tomonaga-Liittinger model should be included. 
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It is also very interesting to compare the results in our paper with similar ones for 
graphene-based nanostructures, where IDEG is provided by chiral-state electrons. A de- 
tailed calculation using RPA for plasmon excitations in both armchair and zigzag edged 
graphene-based nanoribbons was carried out by Brey and Fertig [llj, and only the armchair 
type was found to exhibit undamped plasmon excitations. The most interesting finding in 
Brey and Fertig's work is the presence of metallic nanoribbons with the lowest electronic 
bands given by E±{k) = ±hvF\k\. In these metallic nanoribbons, the electron wave func- 
tions favor forward scattering. The scattering between the branches of different helicity is 
prohibited and corresponding structure factor vanishes. Both helical branches (given value 
of the pseudo-spin) contribute equally to the response . As a result, the polarization re- 
tains time-reversal symmetry. For metallic graphene nanoribbons, we obtain the interband 
polarization function 



(18) 



1 



hcu — hvpq + iO+ hcu + hvpq + iO+ 
2 hvpq'^ 



vr {huj + iO+y - {hvpqY ' 
Note that the polarization subindex here is related to the valley index rather than to the 
pseudo-spin. At the same time, the intraband polarization functions, n_ _(g, u + iO^) and 
n+^+(g, CJ + zO"^), vanish due to time-reversal symmetry. This, in turn, results in the plasmon 



dispersion in Eq. ( 16 ) with electron Fermi velocity independent of the electron density. In 



other words, we have vp = const as in the case of 2DTI. Except for the spin factor of 2, the 



polarization function in Eq. (18) is identical to the two-component HL in 2DTI (i.e., both 



electron spin states are excited). From Eqs. (13), (16) and Fig.|4| one can easily understand 
that the energies of the collective excitations for the two and single component HLs are well 
separated in spectrum, but they share the same particle-hole excitation region. 

For graphene, the plasmon excitations in doped semiconductor armchair nanoribbons are 
similar to the conventional IDEG plasmon dispersion. The localization of electron wave 
functions along two ribbon edges can be obtained in zigzag nanoribbons. However, the 
cone-like band structure disappears for this case. Moreover, the particle-hole excitations 
cover a broad region, instead of a narrow line, and the plasmon dispersion which falls into 



13 



this broad region becomes Landau damped. 



IV. SINGLE COMPONENT HL RESPONSE AS A SPIN-DENSITY WAVE. 

The above paragraph was arguing on strong plasmon branch separation in single and two 
component electron HL. However, single component HL is somehow an artificial concept 
since the topological states always appear as Kramers pair. Below we shall demonstrate 
that the density-density plasmon excitations of single component HL are actually equivalent 
to the spin-density plasmons in two component HL. As a by-product, we shall also consider 
spin-spin waves and compare them to those in a conventional 2DEG. 

Now let us provide theoretical description of the spin-spin and spin-density waves. For 
that purpose, we shall need the inter-spin polarizations (H-t zp) as well as intra-spin polariza- 



tions (H-I- -I-). The intra-spin polarizations are given in Eq. (18). The interspin polarizations 
are also connected by H+__(g,(X') = H_^+(g, w). Here, 

' J E+{k) - E_{k + q) + u ^ ' 

= --^ [Log {2Ak^ + Aq + u) + Log {-2Ak^ + Aq + u)] , 

where we have introduced a cut-off wave vector kc > q due to logarithmic divergence of the 
above integral when the chemical potential is set to zero. In the long wave approximation 
{Aq + uj)/2Akc ^ 1, we can simplify the inter-spin polarization as: 



n-,+(g,c^) (20) 

1 

'2A 



m + Log {2Akc) - {u + Aqf 



Therefore, two-component helical liquid (Kramers pair) obeys the same response function 
as multi-level electron plasma p8] . The generalized nonlocal, dynamic dielectric-function 
matrix connects external Vext and induced Vind perturbation potential matrix elements: 



Ve.tW') = J2C'■,^u'{^^\V^nd\^^') , (21) 



= V,. - V,S;'P,,^,U,,,iq, u) . (22) 
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Here, the superscript p, p indicates density-density response; Vq = V{q) in Eq. (10); the 
sates |p) = |±) indicate the spin component of the Kramers pair (eigenfunctions of Pauh 
matrix). The spin overlap function is given by: 



(23) 



with (To being an identity operator. Exphcitly, Eq. (22) can be written as: 



cP'P 



( 1 - K?n+,+ -i/gn+,+ \ 

10 

1 

V -^n-.- 1 - K,n_,_/ 



(24) 



whose inverse is: 



(25) 



i-i^,(n+,+ + n.,_) • 

In order to consider spin-spin response, we shall introduce the spin-lowering and spin- 
raising operators by their action on the spin states: 



^J+) = |-); ^il-) 
^tl+) = |o); ^tl-) 

(0|±)=0. 



|o) 
1+) 



The spin-spin generalized functions are given by Eq. (22) with 



The explicit form of the dielectric function is: 

e^'-^ = diag[l,l-\/gn+,_,l,l] 
6^'^ = diag[l,l,l-V;n_,+ ,l] 

and the spin-spin waves are given by real zeros of : 



(26) 



(27) 



-1 



1 - V, (n+,_ + n_,H 



(28) 
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FIG. 7: Spin-Spin wave (Eq. (|30[)) and the plasmons (Eq. (25)) in TI. Thick blue curve correspond 



to spin- wave with large critical cut-off 1 / {2A)VqLog(2Akc) ^ 1. Thick green curve is the spin- wave 
with the cut-off in the linear regime of the dispersion 2Akc = 1- Dotted curve is the plasmon mode 
while thin black line is the particle-hole excitations. 

Dispersion curves of the spin-spin waves are shown in Fig.[7]Note that in the limit of large 
cut-off wave vector 1 / {2A)VqLog{2Akc) ^ 1, the spin- wave is analogues to Q- inter-spin 
branch in case of Rashba spin-orbit [18] . The plasmon branch is similar to Qq of intra-spin 
plasmon. n„ inter-spin branch has no analogue in TI case. The reason to that is four 
possible inter-spin transitions in the Rashba split 2DEG as opposite to two of those in TI. 



The spin-density response is defined by Lozovik [7] as Eq. (22) with the spin overlap given 



by: 



(29) 



where a = [ax = crt + (^Ij = (c"t ~ '^i) ) (^z] ■ One can show that 



(30) 



x,x\ — l 



1 - Vq (n+,_ + n_,H 
1 



^ ' i-Vq (n+,_ + n_,+) • 

Due to the fact that in 3DTI the spin wave function component depends on wave vector, 
we are left with only transverse (a^) component of the response survive. In our case, those 
actually provide the spin-spin waves. Those spins may be coupled to the density via e'''^, e'''^. 
However, the inverse of those matrices does not provide any resonances, and then, such 
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coupling may be neglected. On the other hand, we find the rest inverse as: 



1 - V, (n_,_ - Hh 



(31) 



(e^-^)-! \ . (32) 

In the above two equations, the time reversal symmetry is broken, similar (but not identical) 
to the case of single component HL discussed in the previous section. There are four spin- 
plasmon modes: 



7r<±^ = Tqv, ± y^(vrAg)2 + {qV,Y . (33) 

(34) 

Apart from a factor of two, in the long wave approximation we shall obtain the mode given 
by the single component HL: 

= = -2u;ogln(g) . (35) 

TT 



Due to Eq. (30) the combined spin-density response can be written via the spin-density 
overlap as: 

S^^^vv' = if^Wolfi') {i^Wo ± a,\u') , (36) 
so that we recover the single-component HL response: 

{^"''^T' - ■ (37) 

To conclude this paragraph we have demonstrated that the spin-density response mimics 
that of single component helical liquid. Also we found similar behavior of spin-spin waves 
in TI to those provided by Rashba spin-orbit split in a conventional 2DEG. 



V. CONCLUDING REMARKS 



The dispersion of intraband plasmon excitations in a semi-infinite inverted HgTe/CdTe 
quantum well has been derived within the random-phase approximation based on a cal- 
culated edge-locahzed topological state of electrons in a single-component helical liquid 
using the BHZ model. Under the perturbation from a linearly-polarized incident light, the 

17 



unique properties in the collective excitation of these edge-bound electrons with a broken 
time-reversal symmetry has been explored. Our calculations predict the plasmon dispersion 
u)p{q) ~ —ujoQ In(ga) for such a single-component helical state in the long-wavelength limit, 
in sharp contrast with UJp{q) ~ —uiQ^y —ln{qW) found for a one-dimensional electron gas 
in a quantum-wire system. Moreover, uq in our plasmon dispersion is independent of the 
linear electron density, similar to the case for a metallic armchair graphene nanoribbon. On 
the other hand, the plasmon dispersion of the two-component helical liquid is found to be 
identical to that of a armchair graphene nanoribbon except for the spin perfector and a 
characteristic-width scaling of the wave number. The particle-hole excitation region shrinks 
into a straight line in our system, in comparison with a wide region for a conventional one- 
dimensional electron gas. The plasmon energy of the single-component helical state is well 
separated from that of the two-component helical state but they share the common particle- 
hole excitation region in the excitation spectrum. The density-density plasmon excitations 
of single component HL are equivalent to the spin-density plasmons in two component HL. 
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